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Abstract 



The topological pressure introduced by Ruelle and similar quantities describe 
dynamical multifractal properties of dynamical systems. These are important 
characteristics of mesoscopic systems in the classical regime. Original defini- 
tion of these quantities are based on the symbolic description of the dynamics. 
It is hard or impossible to find symbolic description and generating partition 
to a general dynamical system, therefore these quantities are often not ac- 
cessible for further studies. Here we present a new method by which the 
symbolic description can be omitted. We apply the method for a mixing and 
an intermittent system. 
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In recent years the application of the thermodynamic formahsm [|| (TF) in analyzing 
dynamical multifractal properties has been widely accepted. Besides giving an il- 

luminating analogy with the statistical mechanics, it provides a deeper understanding of 
nonanalytic behavior in the scaling properties of trajectories in dynamical systems which 
can be interpreted as phase transitions 0. The topological pressure (TP) plays a central 
role in the TF. The computation of the TP is difficult in general, since the knowledge of 
the symbolic dynamics of the system and its generating partition is inevitable. There 
is no general theory at present, which can provide a systematic and numerically realizable 
method to construct a generating partition to a general dynamical system [|7[. Therefore 
the TP has been computed mostly for low dimensional maps and billiard systems, where 
the symbolic dynamics is accessible, or has been computed with averaging the generalized 
Lyapunov exponents In mesoscopic devices like dots, antidots and wells Hamiltonian 
dynamical systems with smooth potentials play significant role. It would be very enlight- 
ing to apply TF for such non-trivial systems. In this Letter we introduce a new technique 
based on a suitably defined correlation function to measure the TP. The advantage of 
the method is that the detailed knowledge of the system is not required and the calculation 
can be made with low computational demand for either hyperbolic or mixing systems. The 
new method makes possible to measure TP by following one single trajectory of an ergodic 
system. 

The TP for maps can be defined as the logarithm P{q) = logzQ{q) of the leading zero 
zo{q) of the Ruelle zeta function l/({z,q) = Hp (l ~ ' ^^ere the infinite product 
goes for the non-repeating or primitive periodic orbits of the system, Ap is the largest 
eigenvalue of its stability matrix Jp and Up is its period. This definition has been extended 
for continuous fiows in Refs. |]l,|10|l- The TP is the leading zero P{q) = so{q) of the zeta 



function l/({s,q) = Hp (l ~ |T7?) ' where Tp is the time period of the periodic orbit. The 
TP is equivalent with the free energy [|lT| qF{q) and also with the Renyi entropies in 
hyperbolic systems which are defined in a slightly different way in the literature. In the 
definitions above, the problem to find the symbolic dynamics is hidden, since it is necessary 
to locate the periodic orbits in the system. At certain values of q the TP has a special 
meaning [0. For q = the quantity —P{0) is the topological entropy per unit time. For 
g = 1 it is zero for bound systems, it is the escape rate for open systems and its derivative 
gives the metric entropy per unit time. The parameter q where the TP is zero {P{q) = 0) 
yields the fractional part of the fractal dimension of the repeller in scattering systems. The 
time — l/P(l/2) gives a lower bound for the quantum resonance lifetimes in semiclassical 



approximation for open systems ||T3|-0 



Zeta functions can be related to transfer operators. In Ref. ^ we have introduced a 
transfer operator in arbitrary finite dimensions, whose largest eigenvalue is related to the 
TP. Here we briefiy repeat the main idea. 

The dynamical system has to be extended to the tangent space of the fiow, where stability 
of the orbits is multiplicative. We adjoin the (i-dimensional transverse tangent space ^ G 
TUx, ^{x) ■ v{x) = 0, to the ((i+l)-dimensional dynamical evolution space x E U C. R*^"*"^. 
The dynamics in the {x, ^) G t/ x TUrc. space is governed by the system of equations of 
variations |T7| : 



x = v{x), ^ = Dv(x)^. 
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Here Dv(x) is the transverse derivative matrix of the flow. We write the solution as 

xit) = f\xo), ^it) = 3\xo)-^o, (1) 

with the tangent space vector ^ transported by the transverse stability matrix J*(xo) = 
dx{t) / Oxq. In order to determine the length of the vector ^ we introduce a signed norm. An 
example is the function 



g : TU^ ^ R, g 



6- 



(2) 



Any vector ^ G T[4 can now be represented by the product ^ = Au, where u is a unit vector 
in the signed norm: g{u) = 1, and the factor 

A*(xo, uo) = g{m) = 9iJ\xo) ■ uq) (3) 

is the multiplicative "stretching" factor 

A*'+*(xo,uo) = A*'(x(t),u(t)) A*(xo,uo). 

The u evolution constrained to ETg^^, the space of unit tangent vectors transverse to the 
flow V, is given by rescaling of (||): 



u' = RHx, u) 



A* (a;, u 



-J fx) • u. 



(4) 



Eqs. (|T]), (|) and (^) enable us to define a multiplicative evolution operator on the extended 
space U X ETg^^ 



C\x' , u'; X, u) 



6{x' - fix)) 



5{u' - R'{x,u)) 
|A*(x, u)|9-i 



(5) 
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This operator is the generalization of the well-known transfer operator Lq(x, y) =| f'{y) 
6{x — f{y)), for arbitrary dimensions and continuous time. In analogy with the one di- 
mensional problem, we would like to determine its leading eigenvalue in order to get the 
TP. 

Our new method uses the relation of the generalized transfer operators introduced above 
and the correlation decay in chaotic systems. Let's take two arbitrary smooth observables 
A{x, u) and B{x, u) depending on the real variables x and on the tangent space variables u. 
We can define their generalized correlation function as 



C%it) = {A{x,u)\A\x,n) l'-'^ Bif\x),R\x,u))), 
where the average is uniform in the extended space 



{A(x, u)) = J dudxA(x, u). 



(6) 



(7) 
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With the help of the operator the correlation function can be rewritten as 

C%{t) = {A{x, u) J du'dx'C\x, u; x', u')B{x', u')). (8) 

We can introduce the eigenvalues —Sn{q) and eigenfunctions (which are distributions in 
general) ipn{x, u) of the stationary problem of the operator (|^) 

e-*'('^)Vn(a;, J du'dx'C\x, u; x', u')Vn(x', u'). (9) 
For Axiom A systems |]l8l one can expand B{x', u') on the eigenbasis: 

B{x', U ) = bn1pn{x', U). (10) 
n 

The correlation function then becomes 

Clsit) = Ee-^'(«)*6„(A(x,u)V^„(x,u)). (11) 

n 

For large t this sum is dominated by the term corresponding to the largest eigenvalue —So(q) 

Clnit) ~ e-^°(^)*C, (12) 

where the constant is C = bQ{A{x,u.)tpQ{x,u))) . For general, non-Axiom A systems, the 
decomposition ( pUD is not valid, since the eigenfunctions (distributions) not necessarily span 
the full function space. Nevertheless, at long times, the decay is always dominated by the 
largest eigenvalue in accordance with ([T^). 

Therefore the TP can be read off from the semi-logarithmic plot of the correlation func- 
tion for sufficiently large time t. 

For an ergodic system the phase space average (0) can be computed as a suitable time 
average for a sufficiently long ergodic trajectory. If, for the simplicity, we choose the special 
observables A{x,u.) = 1 and B{x,u.) = g{x), where g{x) is the equilibrium probability 
distribution of the system, the correlation function becomes 

cut) = hm i r dt' I A*'+*(xo, uo) (13) 

T^oo i Jq 

where Xq and uq represent an initial condition for a long ergodic trajectory. 

To implement the technique for a two dimensional classical problem, we need a convenient 
choice of the g{C,) function. For 2-dimensional Hamiltonian dynamics let the 2-dimensional 
Poincare section return map be Xj+i = f{xi). The stability matrix of cycle p is a product of 
the 2x2 stability matrices 



J, 



where Aj = dfi{xj)/dxi, and so on. Assume the signed norm (Q) and multiply an initial 
unit vector by the first stability matrix in the product. The resulting vector can be written 

as 
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Hence the dynamics acts on the unit vectors as a rational fraction transformation 

Kk+l = R{Xk, Kk) = — — — , (14) 

UfcKfc + J^k 

with the signed norm of the iterated vector given by 

n 

A"(a;i,Ki) = l[iQn, + Di). 

i=l 

In the case of 2-dimensional bilhards, Kn is the Bunimovich-Sinai curvature pTI]. The 



Bunimovich-Sinai curvature K{t) is the local curvature of the front (horocycle) formed around 
a central trajectory by nearby orbits started from a common point of origin with the same 
energy. For a periodic orbit = kq, the unit vector is an eigenvector of the stability 
matrix, and the corresponding eigenvalue is Ap = Y\{CiKi + Di). The k can be defined also 
for continuous time |2^. Instead of the rational fractional transformation (p!^, we get a 
differential equation for K{t). The "stretching factor" becomes the integral 

A*(a;o, Ko) = exp (^J^ k{kq, xq, t')dt'^ . (15) 

The time evolution of n{t) for 2-dimensional Hamiltonian systems with hamiltonian H = 
\{pi + Py) + U{x,y) can be derived easily: An infinitesimal configuration space volume 
Vq, formed around a central trajectory by nearby orbits started from a common point of 
origin with the same energy, is stretched after time t with the factor A*(xo, hq) orthogonal 

to the central trajectory and with w offlm^ff?^ along the trajectory, due to the change of 



2{E-U{x{0)) 

the velocity along the trajectory. The change of the volume in unit time can be written 
as ^ = X{t)Vt, where the expansion rate is the sum of the orthogonal and the parallel 
stretching rates: 

X{t) = K{t) + j^\\og{2{E-U{xm)- 

Computation of K{t) is then reduced to the problem of finding the expansion rate of volumes. 
This can be recovered by investigating the evolution of trajectories, which deviate infinites- 
imally 6x{t) from the central trajectory x{t). These are described by the linearized Newton 
equation 6x(t) = — D^f/(x(t))(5x(t), where D^U{x{t)) is the second derivative matrix of the 
potential. Since the trajectories are restricted to the same energy surface as the central 
trajectory, we get the additional constraint: x{t) ■ 6x{t) + VU {x(t))6x(t) = 0. The evolution 
in two dimensions can be reformulated to the linear form 

6x{t) = M{t)6x{t), (16) 

where the four elements of the two-by-two matrix M(t) can be constructed from the four 
vector components of 6x{t) and 6x{t). The evolution of the matrix M(t) can be determined 
from the linearized Newton equation: M(t) = —NP{t) — T)^U{x{t)), and the energy conser- 
vation constraint yields: x{t)'M.{t) + VU{x{t)) = 0. The expansion rate of the configuration 
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space volume is the divergence of the velocity field (0), which is the trace of the matrix 
M(t): X{t) =TrM{t). 
For the K,{t) we get: 



. _ _ 2 _ JPydxU - PxdyUf _ Pp^xU - 2p^pyd^yU + vjdyyU 

{Pl+PIY Pl+Pl ^ ^ 



The numerical calculations have confirmed the applicability of the method. For concrete 
calculations we have chosen the Hamiltonian of the Anisotropic Kepler Problem (AKP) 



H = -{pi+pl) 



which describes a Bloch electron with anisotropic mass tensor in a Coulomb potential and 
the Hamiltonian 

h = -{pI+pI) + -xW 

which can be considered as a model of a soft wall mesoscopic cross junction [Q. We have 
measured the correlation function (|13]) for the AKP at parameter a = 5.0 and for the 
potential at energy E = —0.5 up to T = 13000 and E = +0.5 up to T = 30000 respectively, 
for several initial conditions while we integrated the equations with dt = 0.001 using the 



fourth order Runge-Kutta |2J] method with adaptive step-size control. Both systems can 
be considered practically ergodic, since regular islands cover less than 0.005 percent of their 
phase space p5|. The main difference between these systems is that the AKP is a mixing 



one at the chosen value of a, while the other one is intermittent. 

For the AKP the correlation function (|T3|) reached the asymptotic behavior very quickly 
for all q (Fig. lb), while for the intermittent system one can see the effect of the critical 
slowing down at g > 1 (Fig. 2b). The measured TP functions are depicted on Fig. la 
and 2a. 

We hope, that the method presented here can be applied for many djTiamical systems, 
where other methods have failed so far. 

The authors are grateful to P. Cvitanovic, H. H. Rugh, T. Tel and P. Szepfalusy for 
discussions. G. V. thanks the European Community Nonlinear Network the fellowship during 
his stay in Orsay. This work was partially supported by OTKA T17493, F17166 and the 
Foundation for the Hungarian Higher Education and Research. 
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FIG. 1. a) The measured Topological Pressure (per unit time) function for the AKP as a 
function of q with the asymptotes, b) The common logarithm of the measured correlation function 
of the AKP as a function of time t. Different straight lines correspond to different values of q, due 
to the fast convergence to the asymptotic behavior. 
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FIG. 2. a) The measured Topological Pressure (per unit time) function for the ^x^y^ 
as a function of q. b) The common logarithm of the correlation function of the ^x^y^ potential as 
a function of time t. One can observe the critical slowing down of the convergence as q becomes 
greater than 1 
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